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We study electromagnetically induced interference at superconducting qubits. The interaction between qubits 
and electromagnetic fields can provide additional coupling channels to qubit states, leading to quantum interfer- 
ence in a microwave driven qubit. In particular, the interwell relaxation or Rabi oscillation, resulting respectively 
from the multi- or single-mode interaction, can induce effective crossovers. The environment is modeled by a 
multi-mode thermal bath, generating the interwell relaxation. Relaxation induced interference, independent of 
the tunnel coupling, provides deeper understanding to the interaction between the qubits and their environment. 
It also supplies a useful tool to characterize the relaxation strength as well as the characteristic frequency of the 
bath. In addition, we demonstrate the relaxation can generate population inversion in a strongly driving two- 
level system. On the other hand, different from Rabi oscillation, Rabi oscillation induced interference involves 
more complicated and modulated photon exchange thus offers an alternative means to manipulate the qubit, 
with more controllable parameters including the strength and position of the tunnel coupling. It also provides 
a testing ground for exploring nonlinear quantum phenomena and quantum state manipulation, in not only the 
flux qubit but also the systems with no crossover structure, e.g. phase qubits. 

PACS numbers: 85.25.Dq, 03.67.Lx, 32.80.Xx, 42.50.Hz 



I. INTRODUCTION 



As controllable artificial atoms, superconducting Josephson devices provide new insights into the fundamentals of quantum 
mechanics and also hold promise for use in quantum information science as qubits.'"^ Recently, a large number of experiments 
associated with the strong driving have been performed such as Rabi oscillation in the multiphoton regimer^ Landau-Zener- 
Stiickelberg (LZS) interferometry,"^"^ and dressed states in a strongly driven two-state system42i2£ 

In a superconducting flux qubit, when the external flux bias $ « 0.5<i>o, the classical potential energy has a double well 
structure parameterized by S^dc = $ — 0.5$o and the kinetic part of Hamiltonian provides quantum tunnel coupling A between 
two wells. In the millikelvin temperature 10 mK, energy-level quantization emerges and the states localized in each well are 
separated by the energy intervals on the order of the plasma frequency ujp, which is usually much larger than A. The two lowest 
states in each well, |0) and 1 1), form an effective two-level system (TLS), which can be manipulated by application of microwave 
pulses. If the external control parameter is changed so that the qubit traverses the crossover, a non-adiabatic transition between 
two states occurs. This transition is known as Landau-Zener (LZ) transitioni^ When the qubit is subjected to a strong periodic 
microwave field, consecutive LZ transitions between two states at the crossover may result in LZS interference^a2i21 which has 
been realized in two-state systems (TLSs) such as superconducting qubits i ^'-^-'^-'^'^''^ — 

In this article, different from LZS interference, we study the strongly microwave driven qubit in the presence of additional elec- 
tromagnetic interactions based on the general theories of TLS s^^-"'^'-^^" — . The interactions between qubits and electromagnetic 
fields provide additional coupling channels to qubit states, which will detour the evolution path to lead a new interference. We 
name the new interference as electromagnetically induced interference (EII). We choose two typical electromagnetic systems: 
(i) a thermal bath with multi-mode electromagnetic fields and (ii) a weak single-mode field, to discuss interesting phenomena 
and applications. 

In practice, superconducting qubits inevitably couple with electromagnetic environments which can often be described as a 
multi-mode thermal bath of harmonic oscillators,"*" experiencing dephasing and energy relaxation simultaneously. Both pro- 
cesses contribute to decoherence which may hinder qubit performance in solid-state quantum systemsi^^ In general, the bath 
couples to the qubit longitudinally, causing the dephasing related with the environment property at low energies ^ ojp. The 
transverse interaction between the qubit and bath produces the interwell relaxation between two lowest qubit states. Here we 
present a microscopic model, in which the relaxation induced coupling contributes to asymmetrical transitions to form interfer- 
ences, i.e. relaxation induced interference (RII), independent of A under the strong driving. In previous works r^i^ia^ial RII is 
usually overlooked and the interwell relaxation is described empirically, while RII always exists and interacts with the primary 
LZS interference. Moreover, we show that the relaxation can generate population inversion in a strongly driven TLS through 
the interference, which violates the intuitive picture. The model presents deeper understanding to the effect of the environment. 
It also provides a useful tool for characterizing the interaction of qubits and environment, e. g., the relaxation strength and the 
characteristic frequency of the bath. 

Since the interwell relaxation is determined by the environment and hard to manipulate, we turn to a more controllable case. 
Then we consider applying the weak monochrome longitudinal microwave field, which is always used to demonstrate Rabi 
oscillation and serve as a basic method to manipulate quantum states^i^d^^ Rabi oscillation can construct new controllable 
coupling channels which act as effective energy level crossover. As a result, the new interference, which we call Rabi oscillation 
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FIG. 1: (Color online), (a). Schematic energy diagram of a flux qubit. The dotted curve represents the strong driving field A cos cut. The 
field through the tunnel coupling A forms a LZS interference, exchanging photons with the qubit. (b). Quantum tunnel coupling exists 
between states |0) and The interaction between a qubit and an electromagnetic system (such as the environment bath or a single-mode 
electromagnetic field) would form new couplings between the two states. 



induced interference, (ROII) is generated under the strong microwave field. The crossover produced by Rabi oscillation is fully 
tunable from the position to the strength. ROII offers an alternative way to operate qubits in the strongly driven region and the 
potential of further application in the time-dependent interference which can manipulate the qubit state on the Bloch sphere over 
short timescales. In addition, it also provides a testing ground for exploring nonlinear quantum phenomena and quantum state 
manipulation, in the general TLS containing no crossover structure, e. g. phase qubits. 

The paper is organized as follows. In Sec. II, the basic models are introduced first. In Sec. Ill, through the perturbation 
theory, we calculate transition rates based on the relaxation induced coupling. Then we present a discussion about relaxation 
induced interference. In Sec. IV, we discuss the competition of two interferences where one is RII and another is the primary 
LZS interference. By tuning the relaxation strength and the characteristic frequency of environment bath, we show different 
interference patterns in saturated and unsaturated cases. Finally, we compare our model with the phenomenological relaxation 
theory. In Sec. V, we study Rabi oscillation induced interference in the presence of the weak monochrome field. Then we discuss 
new interference mechanism and its application in quantum state manipulation. 



II. BASIC MODEL 



We consider a periodically driven superconducting flux qubit with the two lowest states |0) and |1) in the right and left well 
respectively. (Fig. 1) As we drive the qubit with the microwave <i>ac = ^rf cos tot, where $r/ is the microwave flux amplitude, 
the time-dependent magnetic flux S^{t) = S^dc + ^ac controls the qubit and the time-dependent energy detuning of states |0) 
and |1) can be described as e{t) ~ Eq + Acosut. £o = (l™i| + l^^oD^^dc is the dc energy detuning of states |0) and 
A = {\nio \ + |toi |)(f>r/ is the microwave amplitude in the unit of magnetic flux, and rrii = deo{^)/d^ is the energy-level slope 
of state This qubit can be described by the Hamiltonian^Si^ 

where ax and cr^ are Pauli matrices, az = — |0)(0| and tr-c = |1)(0| + |0)(1|. Firstly we consider the interaction of 

the qubit and the environment bath with a continuum of electromagnetic modes in equilibrium at temperature T (Fig. 1). The 
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Hamiltonian describing the interaction of the qubit and bath reads (we set ft = fcs = 1) 

H = ~—ax —crz-^(Jz y^Tx + Hb, (1) 

where Hb = J^^ii'^l'^i + 1/2) is the bath Hamiltonian, a] and denote the boson creation and annihilation operators 
corresponding to the frequency uji for the bath, Q ^ i'A + a^) is the environment operator, and A; expresses the oscillator 
coupling constant. —tpQax/'i describes the transverse interaction between the qubit and the multi-mode fields (the interwell 
relaxation) that causes the energy exchange between the qubit and bath. —Q(Tz/'2 describes the longitudinal interaction between 
the qubit and the multi-mode fields (the dephasing process) dominated by the low frequency energy from the bath, (p represents 
the transverse interaction strength. Here, we replace the conventional expressions with the noise operator niodela^i21^2i2i^, 
where the environment noise is assumed to follow the Gaussian approximation. As a result all averages can be expressed by the 
spectral density 

1 r°° 

S{u;)^— dte'^ '{Qmm- (2) 
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Now we follow the process in Ref. 33. By taking an interaction picture with respect to Hb, the Hamiltonian is changed to 

e{t) Q{t) A 0g(t) 

Hi = —<7z —<7, - —dx 2~ 

where Q{t) = e*^^*(5e^'^^*. Then we make a transformation to a rotation frame, such that a state vector \'4>) in the interaction 
picture can be expressed as 

IV) =f/o(i)|V"), (4) 

with 

i /•* 

U^{t) = rcxp{-a, [Q{t) + e{T)]dT}, (5) 

where T expresses the time ordering and | V'/) is the state vector in the rotation frame. 
Here we get the effect Hamiltonian 



H'{t) = -l/2cxp{-i I t 
Jo 



e{T)dT} [AU\_ {t)U_ (t) + (t)U\ {t)Q{t)U_ (i)] 1 1) (0| +iJ.c. (6) 

where U± (t) = Texp[±{i/2) Q{T)dT]. Therefore, the time evolution operator in the laboratory frame can be expressed as 

Uiih,t2) - e-'"^*'Uoiti)Uvih,t2)Ulit2)eJ"^*\ (7) 
Uv{ti,t2) = rexp[-i j ' dt'H'it')]. (8) 

Jt2 

Then we obtain the total density matrix of the qubit system plus the bath at the time interval t 

PSBit)^Uiit,0)psB{0)Ul{t,0). 

The system reduced density matrix is p{t) = TiB[psB{t)], where TrB[...] denotes the trace over the environmental degrees of 
freedom. When the system-bath interaction energy is small, we assume psb{0) factors as a direct product psb{0) ~ p{0) ® PB, 
where pB = e^^^^^ is the density matrix of the bath. 

III. RELAXATION INDUCED INTERFERENCE 

Due to the weak relaxation and tunnel coupling compared with the dephasing rate we consider the time interval t within 
which the change of qubit population is slow, but t is much larger than the dephasing time T2 = 1/T2- Then we make the 
approximation with Uv{t, 0) by performing a perturbation expansion. 



f dt'H\t'). (9) 
Jo 



Therefore, the system reduced density matrix can be expressed as 

p{t) = Trs [Ui (<, O)p(O) PS t/J (t , 0)] 

= ^rB[\ I dh I dt2{[At/|(ti)C/-(ii) + 0t/l(ii)g(ti)f/-(ti)]e^'^«'^(")''ni)(O| 
4 Jo Jo 

+H.c.}p{Q) ® pB{[l^Ul{t2) U-{t2) + <i>Ul{t2)Q{t2) U-{t2)]e-'^-" '^^^''^ \1) {0\ + H.c.}]. (10) 



The strong dephasing makes the off-diagonal elements of p{t) decay completely<21 and hence we focus on the diagonal part of 
the density matrix. Then we assume the initial condition p{0) = and obtain poa{t) 

Pooit) = \ I dh f dt2TrBm{[AUlih)U+ih) + (f>ul{h)Qih)U+{h)] 







xe'^»' '^^'>^^\0){1\}\1){1\^ PB{[AUl(t2)U^{t2) + <f)Ulit2)Qit2)U^ 

= \ I dh [ dt2e'^'V^M'^-(A2c/|(i2)C/-(i2)f/l(ti)f/+(ii) + 0'C/|(t2)Q(t2)t/-(t2)f/l(ti)g(ti)t/+(ii) 
Jo Jo 

+0A C/|(t2) U-{t2)Ul{h)Q{h) U+{h) + A0f/|(t2)Q(t2) U-{t2) ul{h)U+{ti)). (11) 
We can calculate the first two average terms in Eq. (1 1) by expanding Uj^{t) to the second order in Q{t), with 

{Ul{t2)U^{t2)Ul{h)U+{h)) « 

1 + i r dt'lT dt {Q{t)Q{t ))- r dt{Q{t)Q{t))] (12) 
Jti Jo Jo 

{ul{t2)Q(t2)U-{t2)ul{ti)Q{ti)U+{ti)) w {Q{t2)Q{ti)) (13) 
Since the last two terms in Eq. (11) are neglected under the average approximation compared with the other two terms, we obtain 

Poo(i) = 4^/ dh f dt2{uUt2)U^{t2)u!_{h)U+(h)) 
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Xe'/'2 + ^ / dtj dt2e'^'2<-^''^Q{t2)Q(h)). (14) 

4 Jo Jo 

The first term in Eq. (14) has been discussed in Ref. 27, where 

{Ul{t2) U-{t2) Ulih) U+{h)) ^ cxp{ J dJs{J) 

' (tl-*2)_i_2jsin "'"r'"' cos iiLlil+M 



X ■ 



-}. (15) 



If the noise is frequency independent (white) in low frequency region tu /27r < 1/t, S{uj ) w 5(0), Eq. (15) can be integrated 
to yield 

([/|(t2);7-(i2)C/l(ii)t/+(ti)) = e-r^l*^-*il, (16) 

where T2 = ttS{0). Then we focus on the second term in Eq. (14) which is related to the interwell relaxation. With the Bessel 
functions 



^i{A/uj) sinujt 



^ J„(A/L.)e™-*, (17) 



where Jn{A/uj) are Bessel functions of the first kind, for the high frequency driving w ^ AJ„(y4/cj), where n is the closest 
integer to eq / uj, we apply the rotating-wave approximation to make 



g2— sinLjti—i— sinutQ 



^j2^^yrMt,-t,)^ (18) 
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where we neglect the nonresonant terms with exponent oscillating quickly with high frequency comparing to the time scale of 
the dynamics t ^ l/[AJ„(A/a;)]. Hence we have 
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■ f dh f dt^Y.Jni^y^''''"'^''^''-''''{Q{t2)Q{h)). (19) 

Jo Jo 
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With the inverse Fourier transformation 

{Q(t2)Q(t{)) = / dc^'e-*'^'(*^-*i)5(w'), (20) 



we find 

a rt 



Poo 



4 Jo Jo 

dh dt2 j dLj'^ j2(^)e*(^°+"'^+"')(*i-*^)5(w'). (21) 



4 Jo 



For the time interval t, which is much larger than the dephasing time, we extend the integration limits t to oo and obtain 

A2 roo ^ poo r , I 

J-oc „ J-oo „ J 

2 (en ~ nuj)'^ + I i 4 w 

where r = ti — <2- Strong dephasing eliminates the off-diagonal elements of quickly and changes the density matrix 
equation to the rate equation. Considering the small time interval, i.e. poQ{t)t 0, we have the rate equation 

Poo(i) = VKio + Tio, (23) 
where quantum tunneling A resulted Landau-Zener transition rate from state 1 1) to |0) is 

and Woi (eo) = Wio{—£o). For the strong low frequency noise with l/w' spectral in low frequency part cj' <C l/r in Eq. (15)^, 
we can derive quantum noise resulted LZ transition rate from state |1) to |0) in a similar way 

= M E ^^(^) exp{-l^^^±^^|±^}, (25) 

n 2 



and Woi(£o) = Wio{—eo), where Sp{t) ~ J dw' ^ (1 — cosw'i), and T2 = y J dio' S{lo'). It should be emphasized that 
LZS interference has been discussed comprehensively in Ref. 27, as a result in order to focus on the physical picture of EII we 
will use Eq. (24) when LZ transition in Sec. Ill and IV is mentioned. On the other hand, seen from Eq. (23), the transition rate 
induced by the interwell relaxation from state |1) to |0) is 

Tio = ^ E Jli-)S{~£o - nuj), (26) 

n 

and roi(eo) = rio(— Eq). It is clear that the dephasing would not influence the transition in Eq. (26). In order to clarify the 
physical picture, we consider an important class of spectral densities^^dli 

1 — e~" /-' 
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The spectral density at negative (positive) frequencies corresponds to the emission (absorption) of photons to (from) the bath. 
Generally, the spectral density describes the ability of energy exchange between the qubit and bath. 

First of all we consider a special situation in Eq. (26), where the microwave amplitude A ^ 0, then Eq. (26) is changed to the 
downward and upward relaxation rates 

Toi = ^^(eo), Tio = ^Si-so). (28) 
Seen from the form of S{uj ), the bath satisfies the Einstein relation with a larger downward relaxation rate 

Tio/roi = e-'"/^. (29) 

roi(rio) corresponds to the process that qubit is spontaneously relaxed down (up) from state |0) to |1) (from state |1) to |0)) 
accompanied by emitting (absorbing) photons with energy Eq to (from) environment bath. Eq. (27) demonstrates that for the low 
temperature T < Q.5ujc there is a maximum downward relaxation rate near the detuning Eq = ±ojc- Then at the characteristic 
frequency lUc, the bath can absorb photons from the qubit fastest. For the high temperature, the downward relaxation rate will 
decrease monotonically with the detuning. From Eqs. (27) and (29), we can know with the detuning approaching to zero, 
relaxation rates become constant and equal while with the detuning larger than Wc, transition rates gradually become zero. 
Therefore even if the temperature does not satisfy the relation T < 0.5wc, there still exist the downward relaxation and smaller 
upward relaxation, mostly near the appropriate detuning Eq = ±Wc- Hence the interwell relaxations can asymmetrically couple 
two states in the qubit at the detuning Eq « ±ojc- Serving as beam splitters, the relaxations at the detuning Eq « — Wc can split 
state |1) with a larger possibility than |0) as indicated in Fig. 2 (a). 

In the presence of microwave, i.e., A 0, the stimulated transition will change the population distribution. Substituting S{oj ) 
in Eq. (27) into Eq. (26), we get 

^'a^ ^2,-4, (£o-na;)e-l^o-""l/"- 
~ — -^"^^^ 1 _ e-(-c-n.)/T ■ (30) 

n 

The term roi(rio) can be understood from Fig. 2, where the qubit state decaying (exciting) processes with roi(rio) are 
shown. For the decaying process in Fig. 2 (a), the qubit emits a photon with energy Eq. Part of this energy is absorbed by the 
microwave field and the residue is absorbed by the bath. If the energy absorbed by the field is larger than £o, i-e. nuj > Eq, 
the bath would emit photons with energy nuj — Eq to the field (not notable in Fig. 2). Fqi reaches maximum when the energy 
Eq is larger than the energy absorbed the field, and the residue energy absorbed by the bath is just Wc. It is similar to the 
exciting process in Fig. 2 (b). Fio reaches maximum when the energy provided by the field is larger than the energy £o, and the 
residue energy absorbed by bath is just Uc- Therefore, the presence of microwave will make the energy exchange process more 
complicate. 

Furthermore, Eqs. (26) and (30) show that the microwave amplitude can modulate the transition through the asymmetric 
couplings to form interference. In order to reveal the physical picture, we further simplify the spectral density in Eq. (27) as 

0,LU ^ UJc- 

It is clear that here the bath can only absorb photons with the characteristic frequency lJc- Moreover there is only downward 
relaxation, which means the splitters at Eq = uJc or — Wc will work only for |0} or |1}, respectively. Substituting Eq. (31) into Eq. 
(26), we obtain transition rates 



01 



4 \ 0, eo / '^c + nijj 

^^^^c / 4{^),Eo = -LUc + nu 
0, eo i= — + nuj. 

At the detuning = i Wc, transition rates with the probability J,^ ( ) are modulated by the amplitude. As shown in 
Fig. 2 (a), when the system is initialized in state the microwave sweeps the state along |1) until it reaches the detuning A 
where = at time ti. Then state |1) is split into the superposition of states |0) and collecting different phases in two 
paths. After evolving independently, they are driven back to A (time t2)- 1 1) is split into |0) and 1 1) again, while |0) is not split. 
The relative phase difference accumulated by states |0) and |1) (marked with the blue region I plus yellow region II) results in 
the interference. We call this case relaxation induced interference (RII). 
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FIG. 2: (Color online) Schematic energy diagram of relaxation induced interference, (a) refers to the transition from state 1 1) to |0) through A. 
(b) refers to the transition from state |0) to |1) through B. The dashed red upward arrows mark the multiphoton absorption from the driving 
microwave field while the dashed green downward arrows mark the multiphoton release to the microwave. The dotted arrows describe the 
energy released or absorbed by the bath. Effectively, the resonant conditions are nuj = e + iOc (from state |1) to |0)) and ncj = e — ljc (from 
state |0) to The blue region I describes effective phase difference contributing to the interference at A. The yellow region II expresses the 
phase eliminated by the bath. The green region III expresses the effective phase difference contributing to the interference at B. 



However, according to the simplified spectral densities in Eq. (31), the second term in Eq. (19) can be changed to 

^ dh dt2 j du j2(^)e'(-o+"-+-c)(ti-t2)^ (33) 

which describes that the multi-mode fields transversely produces a phase accumulation iuJct in the transition from state |1) to 
0). Then the effective phase difference contributing to the interference at A is Aflga (blue region I) since the bath eliminates 
part of the phase (yellow region II). At A, the interference phase lS.Oqa ~ //^ dt{£\Q) — £|i) + Wc), where = —e{t)/2 and 
^|o) = £(i)/2, is modulated by the amplitude. The interference condition A9qa = 27rn or A9qa — ^Trn + tt would make the 
interference effect maximum or minimum. The total phase gained over one period referring to A is 9a ~ / (it(e|o) — e|i) +i^c) ~ 
27r(eo +Wc)/a;. Then the qubit state is driven to the detuning B where £o = and a similar procedure happens with an effective 
phase difference contributing to the interference at B (green region III). The total phase gained over one period referring to B 
is 9b ^ § dt{e\i'j + ujc — £jo)) = 27r(— £o + <^c)/'^- The stable interference can be constructed with the resonant condition 9a 
or: 9b ^ 2TTn. It should be mentioned that this resonant condition is different from the interference condition, playing the role 
of stabilizing the interference effect. Therefore, the resonant transition from state |1) to |0) (from state |0) to |1)) would occur 
with the relation shown in Eq. (32). For the actual spectral density as shown in Eq. (27), there would be additional interference 
resonant peak width which is on the order of Uc- Therefore, in order to form clear interference fringes, the microwave frequency 
should be larger than the characteristic frequency. 

Generally, the interwell relaxation provides additional couplings between states |0) and 1 1) to form extra interferences. For the 
flux qubit, together with the primary LZS interference, the two ones would simultaneously determine the final qubit population. 
In addition, due to asymmetrical upward and downward relaxation transitions, there are interesting phenomena which will be 
discussed in the next section. 
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FIG. 3: (Color online). The stationary population of relaxation induced interference. The pattern is obtained from Eq. (34). (a). The 
characteristic frequency LJc/27r = 0.05 GHz with the temperature 20 mK. Features of population inversion and periodical modulation are 
notable, (b). The characteristic frequency ljc/Stt = 6 GHz with the temperature 20 mK. (c). The characteristic frequency uic/'i.Ti = 0.05 GHz 
with the temperature 2x 10^ mK. In above figures, the driving frequency uj/2-k = 0.6 GHz. 



IV. COMPETITION OF TWO INTERFERENCES 



In order to study RII clearly, firstly we assume there is no crossover structure, i.e. A = 0, to separate this interference from 
the primary one. Then the stationary population in state |0) in RII is 

rio(eo) + roi(£o)' 

In the case a; ^ lOc, seen from Eq. (30), at the detuning Eq = nuj — ujc, the multi-photon resonant fringes of the transition 
rate rio(£o) are observable in Fig. 3 (a). Due to the asymmetric transition rates, population inversion also emerges. The 
modulation of the amplitude results that the interference fringes is the Bessel dependence on the driving amplitude. Moreover, 
the peak width of spectral density yields different fringes from those of the primary LZS interference. In the other case uj < ujc, 
the resonant peak width in RII is larger than ui. Therefore, the individual resonances and Bessel function dependence are no 
longer distinguishable and merge into a continuous band [Fig. 3 (b)]. The population inversion also vanishes. In Fig. 3 (a) the 
temperature is larger than the characteristic frequency while in Fig. 3 (c) the temperature is much smaller than the characteristic 
frequency. Although in the high temperature there is not maximum downward relaxation, it is obvious that patterns in Fig. 3 (a) 
and (c) are similar, which agrees the discussion in Sec. III. 



A. Population dynamics 

Then we return to the dynamics of the qubit interacting with the environment bath. In order to describe the time evolution of 
the qubit population, we employ rate equations, in which the population obey 

-{Woi + ToOpooW + {Wio + rio)/9ii(i) = poo{t) 

Poait) + Piiit) ^ I, (35) 

where Wqi and Wio are defined by Eq. (24), Fqi and Fio are defined by Eq. (30). Therefore, RII and the primary one would 
contribute to the final qubit population together. The population in state |0) is 

M^io(eo) + rio(eo) 

Poo[t) = 



rio(eo) + 2VFio(eo)+roi(eo) 



ftanhf— ) Vt^io(£o) + Tip . ir,o(eo)+2W,o(eo)+ro,(eo)]t .ggx 

^ rio(£o) + 2M^io(£o)+roi(eo)^ 

Here we assume at the initial time the qubit is in equilibrium with its environment. 

If the system dynamics time is long enough that the population is stationary, the population in state |0} can be described by 

^ M^io(£o)+rio(£o) 

rio(£o) + 2M/io(£o) + roi(£o) ■ ^ ' 

RII, dependent on the characteristic frequency Wf,, would greatly influence the final population, and hence the characteristic 
frequency ujc would determine the final population distribution. Fio (roi) is on the order of ojc^^a, which is always smaller 
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FIG. 4: (Color online). Calculated final qubit population versus energy detuning and microwave amplitude, (a). The stationary interference 
pattern in the weak relaxation situation. The parameters we used are the driving frequency a;/ 2 tt = 0.6 GHz, the dephasing rate V2/2-n = 0.06 
GHz, the couple tunneling A/27r — 0.013 GHz, (fy^a — 0.0002 GHz, the temperature is 20 mK, and the characteristic frequency Ucj^T^ = 0.05 
GHz. The periodical patterns of RII can be seen, although not clear, (b). The stationary interference pattern in the strong relaxation situation 
with a<^2 = 0.02 GHz and a;c/27r = 0.05 GHz. Since the relaxation strength is stronger, the periodical interference patterns are more notable, 
(c). The stationary interference pattern in the weak relaxation situation with a.^~ — 0.000002 GHz and ljc/'^tt = 6 GHz. (d). The stationary 
interference pattern in the strong relaxation situation with aip^ — 0.0002 GHz and ljc/2tv = 6 GHz. (e). The unsaturated interference pattern 
in the weak relaxation situation. The system dynamics time t — 0.5 /is. The characteristic frequency ujc/2it = 0.05 GHz. a(f>^ — 0.0002 
GHz. (f). The unsaturated interference pattern in the weak relaxation situation. The system dynamics time t = 0.5 /xs. The characteristic 
frequency LJc/27r — 6 GHz. = 0.000002 GHz. (g). The unsaturated interference pattern in the strong relaxation situation. The system 
dynamics time t — 0.5 fis. The characteristic frequency ujc/2-jt = 0.05 GHz. a<j>'^ = 0.02 GHz. (h). The unsaturated interference pattern in 
the strong relaxation situation. The dynamics time t = 0.5 ^s. The characteristic frequency ujc/2tv = 6 GHz. a(j>^ = 0.0002 GHz. The other 
parameters used in these Figures are the same with those in Fig. 4 (a). 



than 14^10 (Wqi) at the region A > Eq. Therefore the stationary population is dominated by the symmetric transition rates Wio 
and Wqi thus nearly 0.5. At the region A < Eq, Wiq (Wqi) can be neglected so that RII would dominate the population. 

In the case ui ^ cjc, for the increasing relaxation strength, i.e. from (jy^aojc ^ A^/r2 to (jPauJc ~ A^/r2, at the region 
A > Eq population inversion can be more pronounced [Fig. 4 (a) and (b)] while at the region A < Eq the population would be 
nearly 0.5 as indicated by Fig. 3 (a). In the case lu < lUc, with the relaxation strength (jp'aujc increasing, at the region A < Eq 
the population would be nearly as indicated by Fig. 3 (b) while at the region A > Eq the visibility of RII improves, but with 
no population inversion, as shown in Figs. 4 (c) and (d). 

For the short system dynamics time, i.e. t ^ r2/A^, the population in state |0) is described by Eq. (36). The relaxation 
strength determines the time the bath needs to excite the qubit and become more important in determining the qubit population. 
For the weak relaxation strength, i.e. (fy^aujc ^ A^/r2, the bath has no enough time to excite the qubit. With the characteristic 
frequency increasing, there is little difference between interference patterns, as shown in Figs. 4 (e) and (f). For the strong 
relaxation strength, i.e. (fP^auJc ^ A^/r2, RII emerges, and results in clearer difference between interference patterns for 
various characteristic frequencies [Figs. 4 (g) and (h)]. 

Therefore, from the short time dynamics, we can obtain the information of the relaxation strength. Furthermore, from the long 
time dynamics, the effect of the characteristic frequency emerges. 

B. Phenomenological relaxation theory and discussions 

If we take an average approximation with iA/ixJsmijjti — iA/a; sin wi2 in the second term of Eq. (14), the final expressions 
of transition rates in Eq. (26) will be changed to 

Tio = ^5(-eo), Toi = ^5(eo). (38) 
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FIG. 5: (Color online), (a). The unsaturated interference pattern obtained by phenomenological relaxation theory. The system dynamics time 
is 0.5 /is. roi/27r =0.000008 GHz. (b). The stationary interference pattern obtained by the phenomenological relaxation theory. roi/27r = 
0.000008 GHz. (c). Comparison of the results of two theories. The blue dashed line expresses the population in |0) with the characteristic 
frequency a;c/27r = 6 GHz, the detuning eo/Svr — 8 GHz, cj>^a = 0.000002 GHz and the system dynamics time t = 13 fis, the red 
dotted line expresses the population with the characteristic frequency ujc/2'jt = 0.05 GHz, the detuning eo/27r = 8 GHz, (j>^a = 0.0002 
GHz, and the system dynamics time t = 16 us. The black line uses the stationary result with the phenomenological relaxation theory with 
Ti/2n — 0.000008 GHz, and the detuning eo/2n — 8 GHz. (d). The unsaturated interference pattern obtained by phenomenological 
relaxation theory. The evolution time is 0.5 fis. roi/27r =0.001 GHz. (e). The stationary interference pattern obtained by phenomenological 
relaxation theory. Toi/2tt =0.001 GHz. Other parameters are identical with those of Fig. 4 (a). 



Then assuming Fqi = const, we obtain Fio = Toie^^"^'^ . This phenomenological theory is consistent with the existing 
experimental results . Then we substitute Eq. (38) into Eqs. (37) and (36), obtaining the stationary and unsaturated 
population. 

Now we discuss the valid region of the approximation. When the relaxation is weak, i.e. Fqi -C A^/r2, the qubit population 
in state |0) in the short time dynamics {t « r2/A^) can be calculated as shown in Fig. 5 (a). For the model of RII, the weak 
relaxation strength (jPauJc <C A^/r2 in the short time dynamics (t « r2/A^) would lead to the population in state |0) as shown 
in Figs. 4 (e) and (f). There is almost no difference between two theories in the weak relaxation and short time dynamics. 

For the long driving time, the stationary situation will be reached. Shown in Fig. 5 (b) is the population in state |0) obtained 
from the phenomenological theory. It is noticed that the stationary situation for the phenomenological theory agrees with the 
long time dynamics of RII. From the detail comparison in Fig. 5 (c) we find that in an intermediary time scale the phenomeno- 
logical theory is still a good approximation. However, if the time further increases, the difference between two models will be 
remarkable, as shown in Figs. 4 (a), (c) and 5 (b). 

Then we consider the strong relaxation, i.e. Fqi ^ A^/r2. Although the relaxation is in the same order of magnitude 
with A^/r2, it can still be much smaller than r2 mentioned in Sec. III. Shown in Figs. 5 (d) and (e) are the results of the 
phenomenological theory in the short time dynamics {t « r2/A^) and in the stationary case, respectively. Comparing Fig. 5 
(d) with Figs. 4 (g) and (h), we can find the disagreement of two models is large in the short time dynamics. In the stationary 
case, the disagreement is larger, as shown in Figs. 4 (b), 4 (d) and 5 (e). This indicates that the phenomenological theory is only 
a good approximation of RII in the weak relaxation case. 

On the other hand, the results of the phenomenological theory in Eq. (38) follow the Einstein relation with a larger downward 
relaxation rate. If we use Eqs. (24) and (38) in the rate equations Eq. (35), the larger downward relaxation rate would make more 
population exist at lower state and population inversion would not emerge. This is the result of well-known Einstein spontaneous 
and stimulated emission (ESSE) theory. Here, under the strong periodic driving, Eq. (38) is replaced by Eq. (30). The relaxation 
described by Eq. (30) will form the asymmetric transitions and produce population inversion, breaking the Einstein relation, as 
well as ESSE theory. As shown in Fig. 3 and 4, at some detunings, the upward relaxation is larger than the downward one and 
the population at higher state is larger than 0.5. 

Nevertheless, for the weak driving A uj, the transition rates Eq. (30) would be changed to the forms of Eq. (38) and 
Einstein relation in Eq. (29) holds again. Therefore, in the presence of the environment bath ESSE theory is only justified in the 
weak driving. Since RII is not dependent on the crossover, these conclusions are applicable to a general TLS. 
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FIG. 6: (Color online), (a). Schematic energy diagram of a strongly driven flux qubit interacting with a weak single-mode field. The green 
solid curve represents the weak field, forming effective coupling between states |0) and |1). 



V. RABI OSCILLATION INDUCED INTERFERENCE 



In above sections, we investigate RII which is produced by the interaction of the qubit and electromagnetic bath. The bath can 
be considered as multi-mode electromagnetic fields. However, the characteristic frequency and coupling strength are dominated 
by the environment thus uncontrollable, which restrict the applications of EII. Therefore, in this section we will introduce a 
more controllable situation. For a clear physical picture, we apply the monochrome high-frequency weak field with only the 
longitudinal component, which would interact with the qubit as shown in Fig. 6. This case is similar with circuit quantum elec- 
trodynamics (C-QED) where qubits are strongly coupled to quantum oscillatorsr^^ especially the spectroscopyi^— However 
C-QED is also hard to change the parameters and here the applied monochrome weak field can realize more flexible control as 
shown below. 

In the following discussions, although the effect of the inevitable environment bath still need to be considered, for a clear 
physical picture we will replace RII by the phenomenological relaxation theory discussed in Sec. IV, not treating it as multi- 
mode fields as discussed in former sections. Then the Hamiltonian (1) is changed to 

F=-fa.-^..-^..-f..+ff., (39) 

where e{t) = Acos{u]t) describes the weak field, A and u! are the amplitude and frequency of the weak field. Now we make a 
transformation to an interaction picture with respect to Hb — e{t)az/2. Therefore, the Hamiltonian is rewritten to 

H, = -^a, - J„(A/^)e-"^V+ " y E JniA/C.)e"^''']cT-. (40) 

Due to the weak driving, i.e. A < uj,we only consider the case n = ±1 in the sum of the above expression. Now the Hamiltonian 
is 

H2 = —cTz —<yz [e - e Jcr+ [e - e Jcr_, (41) 

where the last two terms would lead to familiar Rabi oscillation^''. Then we repeat the process from Eq. (2) to Eq. (10). 
Considering the initial condition p{Q) = 1 1) (1 1, we have the population in state |0) 

4 Jo Jo 

»Ps{[f/^(t2)f/-(t2)(e-'"*^ -e*"*^)]a+}|0)}. (42) 
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We accept the rotating-wave approximation with the weak and strong fields, respectively, obtaining 

4 Jo Jo 

4 Jo Jo 

If the noise spectral density is dominated by the white low frequency noise, following the process from Eq. (19) to Eq. (22), 
we have 

where T2 has the same definition with Eq. (16). Considering the small time interval, i.e. poo{t) « Wio, we obtain the transition 
rate from state |1) to |0) 

'° 2 ^ (£0 + + w)2 + r2 ^ 2 ^ (£0 + - w)2 + r2 ' ^ ' 

with Woi(£o) = W^oi(— £o)- In order to ensure the validity of Eq. (45), the frequency of the weak field ui needs to be much 
larger than the driving frequency ui of the strong field. 

In the absence of the strong field, i.e. A = 0, the Hamiltonian in Eq. (39) is used to describe Rabi oscillations at the detunings 
£0 = idi with the Rabi frequency AJi{A/oj). Rabi oscillation replaces the primary tunnel coupling A at the detunings £0 = 
with two couplings at the detuning A' (sq = oj) and detuning B' [eq = —ui), Aa' — Ab' = AJi{A/uj). Unlike the interwell 
relaxation, Rabi oscillation can construct symmetrical couplings between states |1) and |0), serving as a beam splitter working 
for both |1) and |0). If w > the detunings A' and B' are well separated. Similar with the primary coupling, under the strong 
field the couplings at A' and B' can form effective crossovers to produce interferences, as shown in Figs. 7 (a) and (b). We call 
this case Rabi oscillation induced interference (ROII). It is important to notice that the position and strength of couplings can be 
controlled by the frequency and amplitude of the weak field. Moreover, the exponent forms in Eq. (43) are similar with that of 
Eq. (33). That means the weak field also produces a phase accumulation iuit, which would influence the interference phase and 
total phase accumulation, just as the multi-mode fields do. 

Now we focus on ROIL First of all, we discuss the transition from state |1) to |0) [Fig. 7 (a)]. When the initial state is the 
microwave sweeps the state along |1) (black solid arrow). As reaching the detuning B', state |1) is split into the superposition 
of states |1) and |0), collecting different phases in two paths. When reaching the detuning A', state |1) is split into |1) and |0) 
again. Evolving independently, states |1) and |0) accumulate the same relative phase difference corresponding to couplings A' 
and B', marked with the region containing the blue region 1, and the yellow region 11. However, the effective phase difference 
contributing to the interference at A' is the blue region 1 since the weak field eliminates part of the phase (yellow region 11). In 
a similar way, the effective phase difference contributing to the interference at B' is the region containing the blue region 1, the 
yellow region 11, and the green region 111. The total phase 6a' gained over one period referring to A' is 27r(£o +Cj)/uj, while the 
total phase 9b' gained over one period referring to B' is 27r(— £0 + Cj)/u}. The stable interference can be constructed with Oa' 
or 6b' = 2im. Therefore, the resonant transition from state |1) to |0) would occur at £0 = noj — cj or £0 = nixi + uj. As the 
amplitude changes, the effective phase will also be modified, resulting in the interference patterns. 

The transition from state |0) to |1) is similar [Fig. 7 (b)]. The effective phase difference contributing to the interference at A' 
is the blue region 1 while the effective phase difference contributing to the interference at B' is the region containing the blue 
region 1, the yellow region 11, and the green region 111. The resonant transition from state |0) to 1 1) would occur at £0 = no; — uj 
or £0 = nix) + Cj. This result is the same with that of transition from state |0) to This consistency makes that there is no 
population inversion in this case. 

The energy conservation is satisfied for these transitions as shown in Fig. 7 (a) and (b). For the transition from state |1) to |0), 
if the transition is through the coupling at A' , besides n photons with frequency cj absorbed, the total energy nu) = £0 + w is 
larger than the energy spacing £0 so that the excess energy Cj flows into the weak field. If the transition is through the coupling 
at B' , besides n photons with frequency cj absorbed, the weak field provides an extra photon with frequency uj. Similarly, for 
the transition from state |0) to the qubit emits n photons with the weak field absorbing or emitting a photon of frequency uj 
through the couplings at A' or B' respectively. 

In order to describe the evolution of qubit population in the presence of two couplings, we employ the rate equations, in 
which the qubit population obeys Eq. (35). But here, Wio and Wqi are described by Eq. (45), with roi/27r = 0.000008 GHz, 
Tio = Toie^^"^^, and the temperature T = 20 mK. Figure. 7 (c) is the calculated qubit population versus energy detuning 
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FIG. 7: (Color online), (a) and (b). Schematic energy diagram of Rabi oscillation induced interference, (a) describes the transition from 
state |1) to |0). (b) describes the transition from state |0) to |1). (c), (d) and (e). The interference pattern of population in state jO) obtained 
from Eqs. (45), (46), and (47), respectively. The parameters used here are aj/27r = 2 GHz, AjGj = 0.9, Vrnjl-K = 0.000008 GHz and the 
temperature is 20 mK. Other parameters of the qubit are identical with Fig. 4 (a). 



and microwave amplitude. To understand this complicate pattern, firstly we consider the only effect of the coupling at A with 
transition rates between states |1) and |0) 

W,, = W,, = ^ (46) 

Substituting Eq. (46) into the rate equations, we have the contour plot of the qubit population in state |0) as functions of energy 
detuning and microwave amplitude, shown in Fig. 7 (d). 

Then we only consider the effect of coupling at B' with transition rates between states |1) and |0) 

^, . VK„ . E , (47) 

2 ^ (£0 + nbj + ijjy- + 

With Eq. (47), we obtain the population in state |0) as shown in Fig. 7 (e). Having addressed two couplings separately, we 
compare the results with Fig. 7 (c) and find that they match very well. The accordance gives a clear physical picture that 
interferences based on couplings at A! and B' construct the final interference pattern. 

Rabi oscillation, resulting from the interaction between the qubit and the weak field, changes the primary coupling at eo = 
to two couplings which are controllable in the position and coupling strength. Then two couplings result into ROIl so as to form 
a complicate structure shown in Fig. 7 (c). If we further increase the amplitude of the weak field, multi-photon Rabi oscillation 
would emerge and more couplings can be formed. 

ROII offers an alternative means to manipulate the qubit. It can be applied further in the time-dependent interference which 
can coherently manipulate the qubit state on the Bloch sphere with a single crossingi^ or continuous drivin g In addition, if 
the strong periodic driving is further replaced to the linear adiabatic passage, the interference is changed to Stark-chirped rapid 
adiabatic passages^^iffi which could serve as elementary logic gates for quantum computing. The adiabatic passages should be 
very slow compared to Rabi oscillation and fast compared to the coherence time. However, in our model, with a single crossing 
or continuous driving, the interference nature'"'""^ allows faster manipulation than adiabatic passages, which is important for the 
short coherence time. 

For EII, the longitudinal monochrome field forms the channels which depend on the tunnel coupling A, such as Rabi os- 
cillation in Eq. (45). On the contrary, the transverse multi-mode fields form the channels which are independent on A, such 
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as the interwell relaxation. Therefore, if we employ a transverse monochrome field, the qubit-field interaction will emerge in 
off-diagonal terms of the Hamiltonian Eq. (1) and Rabi oscillation^^ can still be produced but independent on A. As a result, 
even in the system without the double well structure, such as a phase qubit, ROII can still be produced. 

In a recent experiment, ^'^ two longitudinal strong monochrome fields were applied simultaneously on the flux qubit, just similar 
with the above discussion. Interesting phenomena different from our results were observed. The reason is that the amplitudes 
of two fields were modified together. There does not exist the stable coupling resulted by Rabi oscillation. Furthermore, the 
frequencies in two fields were comparable. Therefore, the phenomenon predicted here does not emerge in that experiment. In 
another experiments^ two longitudinal weak monochrome fields were applied on the flux qubit. In this case, two fields interacted 
with the qubit and the energy transfer described above emerged. But the small amplitude made the interference we predicted 
invisible. Nevertheless, both experiments demonstrate the feasibility of our model. 

VI. CONCLUSION 

We have proposed and comprehensively studied electromagnetically induced interference in superconducting qubits based on 
the couplings formed by different electromagnetic systems. Especially, we focus on a flux qubit with the double well structure. 

First of all, we consider an electromagnetic bath which is familiar and unavoidable. We show the transverse multi-mode fields 
form the interwell relaxation, which can asymmetrically couple two qubit states at the detuning Eq = ±cjc- In the presence of 
the strong driving, the relaxation induced couplings can generate the interference, i.e. RJI. Considering the rate equations, we 
show that the primary LZS interference resulted by A and RII interact together 

Our model provides a deeper understanding to the effect of the environment and shows the way to determine the relaxation 
strength and characteristic frequency of the environment bath. For the weak relaxation, we note that the phenomenological 
relaxation theory can be a good approximation. But for the strong relaxation, the phenomenological theory cannot be applicable. 
Moreover we show that RII can break the ESSE theory under the strong driving. As a result, population inversion can be 
produced in a strong driven TLS by the relaxation and the ESSE theory can only be justified under the weak driving. 

Furthermore, we consider a longitudinal monochrome weak field, which interacts with the qubit through Rabi oscillation. 
We show Rabi oscillation changes the primary tunnel coupling into two symmetric beam splitters at Sq = ±a). Different from 
RII, ROII here can be controlled in the position and strength of the beam splitters by the amplitude and frequency of the weak 
field. ROII provides an alternative way to manipulate qubits in the strongly driven region, and characterize the qubit as well 
as the external fields. It provides the potential of further application in the time-dependent interference which can coherently 
manipulate the qubit state on the Bloch sphere with a single crossing or continuous driving. The models used here can be 
extended to other systems with double well such as quantum dotsr^ or even without double well such as phase qubits. Relative 
experiments^ have been already performed to test some of the predictions of this work. 

VII. ACKNOWLEDGEMENT 

Thanks to Xueda Wen for useful discussions. This work was supported in part by the State Key Program for Basic Researches 
of China (2011CB922104, 2011CBA00205), the NSFC (10725415, 91021003), and the Natural Science Foundation of Jiangsu 
Province (BK2010012). 



* Electronic address: yuyang@nju.edu.cn 

' Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73, 357 (2001). 

^ M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, (Cambridge University Press, Cambridge, 2000). 
' J. E. Mooij, Science 307, 1210 (2005) 

J. Clarke and R K. Wilhelm, Nature (London) 453, 1031 (2008) 
^ J. Q. You and R Nori, Phys. Today 58, No. 11, 42 (2005) 
^ G. Wendin and V. S. Shumeiko, Low Temp. Phys. 33, 724 (2007). 

M. H. Devoret and J. M. Martinis, Quantum Inf. Proc. 3, 163 (2004). 

Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Phys. Rev. Lett. 87, 246601 (2001) 
' S. Saito, T. Meno, M. Ueda, H. Tanaka, K. Semba, and H. Takayanagi. Phys. Rev. Lett. 96, 107001(2006) 

W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. Levitov, and T. P Orlando. Science 310, 1653(2005) 
" M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and P Hakonen. Phys. Rev. Lett. 96, 187002(2006) 

D. M. Berns, W. D. Oliver, S. O. Valenzuela, A. V. Shytov, K. K. Berggren, L. S. Levitov, and T. P Orlando. Phys. Rev. Lett. 97, 

150520(2006) 

" D. M. Berns, M. S. Rudner, S. O. Valenzuela, K. K. Berggren, W. D. Oliver, L. S. Levitov, and T. P Orlando. Nature (London) 455, 51(2008) 



15 



M. D. LaHaye, J. Suh, P. M. Echtemach, K. C. Schwab, and M. L. Roukes, Nature 459, 960 (2009). 

G. Z. Sun, X. D. Wen, B. Mao, Z. Y. Zhou, Y. Yu, P. H. Wu, S. Y. Han, Nature. Commun. 1:51 doi: 10.1038 / ncommsl050 (2010) 
J. R. Petta, H. Lu, and A. C. Gossard, Science 327, 669 (2010); 
" G. Z. Sun, X. D. Wen, Y. W. Wang, S. H. Cong, J. Chen, L. Kang, W. W. Xu, Y. Yu, S. Y. Han, and P H. Wu. Appl. Phys. Lett 94,102502 
(2009) 

Y. W. Wang, S. H. Cong, X. D. Wen, C. Pan, G. Z. Sun, J. Chen, L. Kang, W. W. Xu, Y. Yu, and P H. Wu, Phys. Rev. B 81, 144505 (2010) 
" C. M. Wilson, T. Duty, F. Persson, M. Sandberg, G. Johansson, and P Delsing, Phys. Rev. Lett. 98, 257003 (2007) 
C. M. Wilson, G. Johansson, T. Duty, F. Persson, M. Sandberg, and P Delsing, Phys. Rev. B 81, 024520 (2010) 

L. D. Landau., Phys. Z. Sowjetunion 2, 46 (1932); G. Zener, Proc. R. Soc. London. Ser. A 137, 696 (1932); E. C. G. Stiickelberg, Helv. 
Phys. Acta 5, 369 (1932) 

A. V. Shytov, D. A. Ivanov, and M. V. Feigel'man, Eur. Phys. J. B 36, 263 (2003) 
S. N. Shevchenko, S. Ashhab, and E Nori, Physics Reports 492, 1 (2010) 

S. N. Shevchenko, S. H. W. van der Ploeg, M. Grajcar, E. Il'ichev, A. N. Omelyanchouk, and H.-G. Meyer, Phys. Rev. B 78, 174527 (2008). 
A. Izmalkov, S. H. W. van der Ploeg, S. N. Shevchenko, M. Grajcar, E. Il'ichev, U. Hiibner, A. N. Omelyanchouk, and H.-G. Meyer, Phys. 
Rev. Lett. 101, 017003 (2008). 

L. J. Du, M. J. Wang, and Y. Yu, Phys. Rev. B 82, 045128 (2010) 
^' L. J. Du, and Y. Yu, Phys. Rev. B 82, 144524 (2010) 

J. Bylander, M. S. Rudner, A. V. Shytov, S. O. Valenzuela, D. M. Bems, K. K. Berggren, L. S. Levitov, and W. D. Oliver, Phys. Rev. B 80, 
220506 (2009) 

^' S. O. Valenzuela, W. D. Oliver, D. M. Berns, K. K. Berggren, L. S. Levitov, T. P Orlando, Science 314, 1589 (2006) 
^° W. D. Oliver, S. O. Valenzuela, Quant. Info. Proc. 8, 261 (2009). 

^' M. S. Rudner, A. V. Shytov, L. S. Levitov, D. M. Berns, W. D. Oliver, S. O. Valenzuela, and T. P Orlando, Phys. Rev. Lett. 101, 190502 
(2008). 

M. H. S. Amin, and D. V. Averin. Phys. Rev. Lett. 100, 197001(2008) 
" M. H. S. Amin and E Brito, Phys. Rev. B 80, 214302(2009). 
I. Goychuk and P Hanggi, Adv. Phys. 54, 525 (2005). 
M. Grifoni and P Hanggi, Phys. Rep. 304, 229 (1998). 

S. Ashhab, J. R. Johansson, A. M. Zagoskin, and E Nori, Phys. Rev. A 75, 063414 (2007). 
^' Y. Gefen, E. Ben-Jacob, and A. O. Caldeira, Phys. Rev. B 36, 2770 (1987). 

S. N. Shevchenko, A. S. Kiyko, A. N. Omelyanchouk, and W. Krech, Low Temp. Phys. 31, 569 (2005). 

K. Saito, M. Wubs, S. Kohler, Y. Kayanuma, and P Hanggi, Phys. Rev. B 75, 214308 (2007). 
^ A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P A. Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 

L. Tian, S. Lloyd, and T. P Orlando, Phys. Rev. B 65, 144516 (2002). 

C. H. VanderWal, EK. Wilhelm, C. LRM. Harmans, andLE. Mooij, Eur. Phys. IB31, 111 (2003). 

R. W. Simmonds, K. M. Lang, D. A. Hite, S. Nam, D. P Pappas, and J. M. Martinis, Phys. Rev. Lett. 93, 077003 (2004). 

K. B. Cooper, M. Steffen, R. McDermott, R. W. Simmonds, S. Oh, D. A. Hite, D. P Pappas, and J. M. Martinis, Phys. Rev. Lett. 93, 180401 

(2004). 

J. M. Martinis, K. B. Cooper, R. McDermott, M. Steffen, M. Ansmann, K. D. Osborn, K. Cicak, S. Oh, D. P. Pappas, R. W. Simmonds, and 

C. C. Yu, Phys. Rev. Lett. 95, 210503 (2005). 

Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 786 (1999). 
**' D. Vion, A. Aassime, A. Cottet, P Joyez, H. Pothier, C. Urbina, D. Esteve, and M. H. Devoret, Science 296, 886 (2002). 

Y. Yu, S. Y. Han, X. Chu, S.-I. Chu, and Z. Wang, Science 296, 889 (2002). 

J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, Phys. Rev. Lett. 89, 117901 (2002). 
^° I. Chiorescu, Y. Nakamura, C. J. P M. Harmans, and J. E. Mooij, Science 299, 1869 (2003). 
^' J. E. Mooij, T. P Orlando, L. Levitov, L. Tian, C. H. van der Wal, and S. Lloyd, Science 285, 1036 (1999). 

J. M. Martinis, S. Nam, J. Aumentado, K. M. Lang, and C. Urbina. Phys. Rev. B 67, 094510 (2003) 
" I. Chiorescu, P. Bertet, K. Semba, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, Nature (London) 431, 159 (2004) 

A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature (London) 

431, 162(2004); 
" R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 (2008). 

A. Fedorov, A. K. Eeofanov, P Macha, P Forn-Dl'az, C. J. P M. Harmans, and J. E. Mooij, Phys. Rev. Lett. 105, 060503 (2010) 
" J. Johansson, S. Saito, T. Meno, H. Nakano, M. Ueda, K. Semba, and H. Takayanagi, Phys. Rev. Lett. 96, 127006 (2006) 

F. Deppe, M. Mariantoni, E. P. Menzel, A. Marx, S. Saito, K. Kakuyanagi, H. Tanaka, T. Meno, K. Semba, H. Takayanagi, E. Solano and 

R. Gross, Nature Physics 4, 686 (2008) 

M. M. T. Loy, Phys. Rev. Lett. 32, 814 (1974). 

L. E Wei, J. R. Johansson, L. X. Cen, S. Ashhab, and E Nori, Phys. Rev. Lett. 100, 1 13601 (2008). 
" G. Z. Sun, X. D. Wen, B. Mao, Y. Yu, J. Chen, W. W. Xu, L. Kang, R H. Wu, and S. Y. Han, arXiv: 1010.5897v3 



